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We discuss the photovoltaic effect at a p-n hetero junction, in which the illuminated side is a doped 
Mott insulator, using the simplest description of a Mott insulator within the Hubbard model. We find 
that the internal quantum efficiency of such a device, if we choose an appropriate narrow-gap Mott 
insulator, can be significantly enhanced due to impact ionization caused by the photoexcited "hot" 
electron-hole pairs. Namely, the photoexcited electron and/or hole can convert its excess energy 
beyond the Mott-Hubbard gap to additional electrical energy by creating multiple electron-hole 
pairs in a time scale which can be shorter than the time characterizing other relaxation processes. 
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I. INTRODUCTION 

Solar cells based on conventional single-gap band- 
semiconductors have low internal quantum efficiency for 
conversion of solar into electrical energji The main 
reason is that the photon energy absorbed by a single 
electron-hole pair with energy above the semiconductor 
band gap is lost in heat through electron-phonon scat- 
tering and phonon emission through which electrons and 
holes relax to their band edges. The solar energy spec- 
trum contains photons with energies in a rather broad 
energy range compared to any particular semiconductor 
energy gap. If instead semiconductors with a relatively 
large band-gap are chosen in order to avoid the above 
mentioned problem, this cuts off most of the energy spec- 
trum of the solar radiation. Using a stack of cascaded 
multiple p-n junctions with various band gaps is an ex- 
pensive approach to build efficient solar cells. 

A variety of other ideas have been suggested in order to 
increase solar cell efficiency^—. One particular approach 
which increases the photo-current is one in which the 
energetic carriers produce additional electron-hole pairs 
through impact ionization^— at a much faster rate than 
the rate characterizing the relaxation process through 
phonon emission. In order to achieve this, it has been 
suggested thal£ using confining geometries, such as quan- 
tum wells, quantum wires, quantum dots, superlattices 
and nanostructures, the relaxational time scales can be 
significantly affected^— and this allows the possibility 
for impact ionization. Ross and Nozik^ have argued that 
an efficiency of approximately 66%, much larger than 
Shockley-Queisser efficiency limiti for photovoltaics, can 
be theoretically achieved using an unusual inversion, in 
which the chemical potential of the excited electronic 
band is below that of the ground band. Other ideas 
to increase solar-cell efficiency by utilizing the possibil- 
ity of carrier multiplication^— via impact ionization by 
removing and isetropically cooling the carriers through 
special membranes were proposed by Wurfeli^. 

In this paper we discuss the photovoltaic effect for ma- 
terials which are based on doped Mott insulators, as op- 



posed to the familiar semiconductors which are band- 
insulators. A variety of transition metal oxides that are 
predicted to be conductors by band theory (because they 
have an odd number of electrons per unit cell) are, in 
fact, insulators. Mott and PeierlaiS, predicted that this 
anomaly can be explained by including interactions be- 
tween electrons. Today, this can be described within the 
simple Hubbard model, as arising from the strong on-site 
Coulomb interaction. In this case the half-filled band 
leads to an insulator; the added holes or added extra 
electrons beyond half filling lead to conduction. In this 
case also we can imagine that we can form a p-n junc- 
tion separating a region of hole-doped from a region of 
electron-doped Mott insulator. 

There are recent reports^— where the photovoltaic 
effect is observed in a p-n heterojunction of a doped mag- 
nanite (a doped Mott-like insulator) and doped strontium 
titanate (a doped band insulator). In addition, there 
are reports of the photoelectric effect observed on het- 
erojunctions of doped magnanites and doped silico n 24 i 25 
and on a heterostructure^ of YBaiCuzC>T- x on Nb 
doped SrTiO?,. Here, we demonstrate that in the case 
of narrow- gap Mott insulators the photovoltaic effect can 
lead to solar cells of high quantum efficiency, where a sin- 
gle solar photon can produce multiple electron-hole pairs. 



II. MODEL AND CONCEPTS 

In the simplest model of a Mott insulator, a single con- 
duction band is formed from one localized orbital (which 
is, in general, a linear combination of atomic orbitals) 
per unit cell (or site). In the case of half-filling, i.e., 
when one electron per site occupies the band made from 
this orbital, the material is a Mott insulator. In Fig.QJa) 
the half-filled state is shown where the motion of an elec- 
tron from one site to the nearest-neighbor site leads to 
an empty site (or hole) and a doubly-occupied site (DO) 
(Fig.QJb)); this state, because of the DO, is characterized 
by a high energy cost, which is denoted by the so-called 
Hubbard on-site Coulomb repulsion energy U , and this 
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effectively prohibits electron motion through the band. 
In order to have a concrete framework to facilitate the 
discussion we consider the so-called Hubbard model. 



H — H t + Hjj , 
H t = -t (c] a c la - + h.c), 

<ij> 

Hu = U 22 n^n^, 



(1) 

(2) 
(3) 



where rii a — c\ a Ci a and c\ a creates an electron of spin a 
at site i. 
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FIG. 1: Representative states from the SJ, subspaces. 

First, let us consider the limit of U 3> t, because in 
this case there is a well-defined concept for the origin 
of the Mott-insulator gap and this simplifies the concep- 
tual part of the discussion. If we assert that the Hjj 
term is the unperturbed part and we treat H t as a per- 
turbation, the eigenvalues of the unperturbed part are 
equally spaced by an energy U. For the case of half- 
filling, i.e., in an TV-site system with N electrons, each 
energy level corresponds to a massive number of states 
forming subspaces Sq characterized by the same eigen- 
value E n = nU, where n is the number of pairs of DO- 
holes (DO-H). Figs. QJa),(b) and (c) illustrate examples 
of states n = 0, 1 and 2 pairs of DO-H. The subspaces 
Sq are spanned by the states obtained by all allowed 
spin rearrangements. In a Mott insulator, the Coulomb 
energy cost for two electrons to occupy the same orbital 
leads to an insulating state at half-filling. Replacing some 
atoms in the lattice with different atoms, which act as 
donors or acceptors of additional electrons beyond half- 
filling thereby creating DO sites (n-type doping) or holes 



(p-type doping), leads to a conducting state with an ef- 
fective carrier density given, respectively, by the density 
of DO or hole sites. In a doped system, i.e., an TV-site 
system with TV + m or N — m electrons, the subspaces 
S^j or 5™ m are spanned by m DO or holes and n DO- 
H pairs. Examples for S 1 ™ and S 1 1 1 (for n = 0, 1,2) are 
illustrated in Figs. [TJd), (e), and (f), and in Figs. [Tig) 
(h), and (i). For the rest of the paper, we will denote by 
U* the renormalized gap, namely, the energy required to 
promote a particle and to create a DO-H pair starting 
from the half-filled lattice and by W the bandwidth of 
the quasiparticle or quasihole dispersion. 

In the half-filled case and at low enough temperature, 
the model exhibits antiferromagnetic order—, which, 
however, will not play a significant role in our discussion, 
because the Mott-insulator phase can exist at tempera- 
tures higher than the Neel ordering temperature. 
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FIG. 2: An interface between a p- and an n-doped Mott insu- 
lator before carrier diffusion across the interface takes place. 
Some of the DO sites can be relieved from the energy cost due 
to Coulomb repulsion via diffusion of the extra electron to the 
p-doped region to fill a hole. This process stops when the elec- 
tric field, which gradually develops across the interface due to 
this charge redistribution, becomes strong enough to oppose 
further carrier diffusion. The idea discussed here is also ap- 
plicable when the illuminated side is a doped Mott-insulator 
and the other side a doped band-insulator. 

The quasiparticle (hole or DO) spectral function and 
dispersion relation, as well as the spectral weights, have 
been computed using the t — J model^ (which is the 
strong coupling limit (U 3> t) of the Hubbard model with 
J = At 2 /U) within the so-called non-crossing approxima- 
tion for 2D. In this paper we have carried out the same 
calculation for 3D and the conclusions are qualitatively 
similar to the 2D case. There is a well-defined quasi- 
particle which corresponds to a low energy peak in the 
spectral function, and we will use the operator which, 
by acting on the undoped state |0), creates this quasipar- 
ticle with energy dispersion and residue Z^. In this 
limit, while the bandwidth of the quasiparticle dispersion 
£"k is of the order of J, the width of the density of states 
is of the order of ltd (where d is the dimensionality). 

Now we wish to consider a p-n hetero junction as il- 
lustrated in Fig. [5] based on an interface between two 
narrow-gap Mott insulators, one of which is p-doped (left 
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side of the interface) and the other which is n-doped 
(right side) . Because of the energy gain that is achieved 
when a hole of the p-doped material is filled by an elec- 
tron from a DO site on the n-doped part, diffusion of elec- 
trons from the n-doped side will occur and this leads to 
a positively charged layer in the n-doped side adjacent to 
the interface and, vice versa, a negatively-charged layer 
in the p-doped side. The potential energy difference acts 
as a unidirectional collector of the DO and holes created 
upon the incident light. 




(a) 
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FIG. 3: See the text for explanation. 



III. NARROW-GAP NARROW-BAND CASE 

Let us, first, consider a narrow-band narrow-gap Mott 
insulator. An incident photon of energy several times the 
Mott-gap U* has enough energy so that several electron- 
hole pairs could be created. In practice, since the range 
of the solar spectrum is from 0.5 — 3.5 eV, we can consider 
a Mott-insulator with a gap of the order of 0.5 — 1.0 eV 
and a bandwidth in the same range. We need a narrow- 
gap Mott insulator in order to increase the efficiency, 
but on the other hand the gap should be large enough 
to avoid decay of the excited quasiparticles (by phonon 
emission, spin-wave emission or through other possible 
excitations) . 

If the energy Tiuj q of the solar photon is between 
U* < Tujj q < U* + W, the incident photon can create 
a single electron-hole pair by promoting an electron to 
the first Mott-Hubbard band. Fig. [3](a) indicates a two- 
step process in which the incident photon can create a 
state with two pairs of DO-H, if its energy is higher than 
2U*. In Fig.[3J(a), within each atom we depict two atomic 



levels, one which forms the Mott-Hubbard band, and 
a higher level which forms a high-energy band. First, 
the incident high-energy photon is absorbed by an elec- 
tron which hops to the high-energy band as indicated in 
Fig. Eta). Subsequently, the charged atom, which carries 
the extra electron, transfers part of the excess energy via 
the Coulomb interaction to an electron of a different site 
which causes it to hop to a n.n. atom, while itself falls 
to the lowest level already occupied by one electron; this 
leads to the creation of two DO-H pairs. Fig. [31(b) il- 
lustrates the same process in momentum space. If the 
energy of the photon is higher than 2U* , the electron 
(or the hole) can be temporarily photo-excited to a high- 
energy band which is different than the Mott band. The 
electron, for example, can be photo-excited to the state 
|k v) of energy E denoted by the label k, v in the first 
term of Fig. [3jb) . This electron can decay quickly into 
two quasiparticles plus one quasihole state of the Mott 
band, provided that the energy of the initial and fi- 
nal state match. This is caused by the matrix element 
V = (k v, h //| V"|p n, k — p + h /i) of the Coulomb inter- 
action which couples the state |k v), where the electron 
was initially excited by the incident high-energy photon, 
and the states |h /i), |p /i) and |k p h /i), which 
form the Mott band /i. Now, we would like to estimate 
the transition rate for an electron or hole of momen- 
tum k from a high-energy band to undergo the decay 
processes illustrated in Fig. [3jb). Namely, Fig. [3j(b) il- 
lustrates the following two processes: (a) the initially 
photo-excited particle can decay into two-quasiparticles 
plus one-quasihole given as a ha^ p al._ p _ h \0) (first term 
in Fig. [3Jb)), or (b) the initially photo-created hole 
state can decay into two-quasiholes plus one-quasiparticle 
given as a^apCtk-p-hlO) (second term in Fig. [3Jb)). For 
simplicity, the Mott-band index fi has been omitted as 
a label of the a's. Using Fermi's golden rule we can ob- 
tain the decay rate of the bare electron (or hole) into 
two-quasiparticles plus one quasihole as: 

r k (e) = ^ l M l 2,5 ( e - ^ k -p+ h - e p + ( 4 ) 

p,h 

where e = E — 2U*, and M stands for the inter- 
action matrix element which couples the initially ex- 
cited particle state c^JO) (|0) is the ground state) to 
which the incoming solar photon couples directly, and 
the final two-quasiparticles and one-quasihole given as 

a ha P ak_p_hl°)> Le -' (°\ a h a p a k- P +hVcl jO). As can 
be inferred from the main contribution to the decay of 
the photo-excited high-energy electron-hole pair into two 
DO-hole pairs is due to the two processes shown by the 
diagrams in Fig. [3jb) . Namely, when the Mott electron 
absorbs most of the photon energy and it is excited to 
a high-energy band, it becomes possible for the electron 
to decay into two DOs and one hole in the Mott band. 
When the created hole absorbs most of the photon en- 
ergy, i.e., it corresponds to an excitation of an electron 
from an inner valence state to the Mott conduction band, 
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it becomes possible for the hole to decay into two holes 
from the Mott band and one DO state from the Mott 
band. Each case can be approximated by the calculation 
of the decay rate of the high-energy electron or the decay 
rate of the high-energy hole because the other particle 
(the hole or electron respectively), in this order in per- 
turbation theory, behaves only as a spectator. The later 
rates are given by the above Eq. 2J The matrix elements 
M in this expression is given by 



V 

M — —ZhZpZk-p+h, 



(5) 



where the Z's are the quasiparticle residues and N is the 
total number of lattice sites. The effect of the Hubbard 
interaction U has been approximately taken into account 
by the Z, where Z p — (0|a p |p n), and by using the renor- 
malized Mott-band dispersion E^. In order to obtain an 
order of magnitude estimate for T, we will take V to be of 
the order of 1 eV. In a Mott insulator, where the atomic 
limit is appropriate, the charged highly-excited atom in 
the intermediate state of Fig.[3ja) is expected to interact 
rather strongly with the electrons of the other atoms in 
its neighborhood. 

This decay rate can be estimated by the following fur- 
ther approximation 



V k (e)~^-V 2 Z 6 D(k,e), 



(6) 



where Z is the average quasiparticle/quasihole spectral 
weight throughout the entire Mott-Hubbard band and 



D(k,e) 



Lj2S(e-E k _ p _ b -E p +E b ), (7) 



N 2 



P ,h 



the two-particle/one- hole density of states. The approx- 
imation of Eq. [6] is obtained by using the average Z 
throughout the integration. 

The two-particle/one-hole density of states can be il- 
lustrated using an effective tight-binding dispersion, i.e., 

Ey = — 2t*X)u=i C0S (^M a ^)' wnere d is the number of 
dimensions and k^, a^, are the component of the mo- 
mentum and the lattice constant along the direction /j,. 
In Fig. |31 we show D(k, e) for a square lattice. Notice, 
that it is non-zero in an energy range of the order of the 
single-particle/single-hole bandwidth W — 8t* and the 
value of the density of states is of the order of 1/W. The 
rate averaged over all momenta in the tight-binding band 
can be written as 



f(e) 
A(e) 



2tt. 



V 2 Z b A(e), 



(8) 
(9) 



which is the density of states given in Fig. [4] for two- 
dimensions (2D) averaged over all momenta k. Note that 
A(e) is of the order of 1/W in the energy range where 
there are available states. 
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FIG. 4: The two-particle/one-hole density of states for the 
simple tight-binding-like dispersion E(k) = — 2t* (cos(fc x a) + 
cos(k y a)) for a square lattice. 



The value of the residue Z has been estimated in the 
strong-coupling limit using the t — J model in Ref. [28| 
for 2D. We have carried out the same calculation in 3D 
and we found that the value of Z is somewhat larger 
than in the 2D case, because the role of quantum fluc- 
tuations is less important in 3D. Because we are inter- 
ested in obtaining an approximate estimate of T, we 
will use the results of Ref. 28 : for 2D. Since the t — J 
model is the strong-coupling limit of the Hubbard model 
with J/t = At/U, we need to choose U ^> t. Choosing 
J/t = 1 (i.e., U = At) from Table IV of Ref.Hwe obtain, 
W = 0.96i and Z = 0.53. Thus, taking V = U = 1 eV 
yields t = 0.25 eV, W = 0.24 eV and r ~ 0.9 x 10 15 sec" 1 
(using A(e) ~ 1/W). As can be inferred from the trend 
in Table IV of Ref. H|| larger values of t/U will give 
larger values both for Z and for the bandwidth W (sim- 
ilar to the value of the gap as needed) and consequently 
a larger T. However, while we believe this trend to be 
qualitatively correct, the regime At/U > 1 is not acces- 
sible by the strong-coupling limit discussed in Ref. [2Sl . 
Nevertheless, the above estimate indicates that this pro- 
cess can be significantly faster than the decay processes 
through phonor*2i£~— or spin-wave excitations which are 
expected to be of the order of 10 12 — 10 13 sec -1 . In addi- 
tion, the Mott-Hubbard gap and W can be tuned to be 
several times smaller than the range of the energy of so- 
lar photons and, thus, a photon has the energy required 
to create more than one quasi-electron-hole pairs. 



IV. NARROW-GAP WIDE-BAND CASE 

Now, let us consider the case of a narrow-gap but a 
wide-band Mott-insulator, i.e., the case where W 3> U* . 
It is believed and supported by calculations that for any 
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value of W we can tune the value of U to provide a small 
value of the gap U* . It is also known that the value of U 
required to produce a small gap should be greater than 
W/2. Therefore, if we choose W = 8 eV and U* = 0.5 eV 
we will need a U greater than 4 eV. The exact value of 
U is not required for the following calculation, suffice it 
to know that it is greater than 4 eV. When the incident 
photon promotes an electron from the lower Hubbard- 
band to the upper Hubbard band, the energy E of the 
excited electron can be large compared to the gap, say 
E > nil*, and, thus, it can decay into n quasiparticlcs 
plus re — l quasi- holes (n = 2,3,..). The main process 
is depicted by the diagram of Fig. |3{b), where all the 
lines in this case denote quasiparticles/quasiholes of the 
same (Mott) band (i.e., v = fi), and the interaction is the 
Hubbard on-site Coulomb interaction U instead of the V 
discussed previously. 

In this case the decay of the high-energy DO-hole pair, 
produced by the solar photon, into two DO-hole pairs 
is due to the two processes shown by the diagrams of 
Fig. EJb) where the interaction is U instead of V and the 
decay rate is given by Eq. 3] with the matrix element M 
given by Eq. [5] where V is replaced by the Hubbard on- 
site Coulomb repulsion U. Namely, (b) when the electron 
absorbs most of the photon energy and it is excited to the 
upper-Hubbard band to become a high-energy DO site, it 
becomes possible for the DO to decay into two DO states 
and one hole, (b) When the created hole absorbs most 
of the photon energy, it becomes possible for the hole 
to decay into two holes in the Mott band and one DO 
state in the Mott band. Here also, either (a) or (b) can 
be estimated by the calculation of the decay rate of the 
high-energy DO or the decay rate of the high-energy hole 
because the other particle (the hole or DO respectively), 
behaves only as a spectator. 

Therefore, an estimate of the decay rate is given by 
Eq. 2] by replacing V with U which can be approximated 
by HT ~ 2ttU 2 Z 6 /W (See discussion leading to Eq. [5] 
starting from Eq.HJ. Using this expression and W = 8 eV 
and U = 4 eV, the lower bound required to yield U* = 0.5 
eV, we also obtain a large T ~ 0.5 x 10 15 sec _1 (taking 
the same value of Z used in the previous paragraph). 



V. ISSUES REGARDING EFFICIENCY 

A. Impact ionization, Auger recombination and 
other processes 

Impact ionization increases the number of charge car- 
riers per absorbed phonon. However, Auger recombina- 
tion which is the inverse process should not be omitted. 
First, let us begin our discussion from the open-circuit 
situation. We found that for the type of Mott insula- 
tors discussed in the previous section, the characteristic 
time scale for thermalization of the photo-excited elec- 
trons and holes is much shorter than the time needed for 
phonons to thermalize the electronic system with the lat- 



tice. Therefore, it is reasonable to assume that the photo- 
created DOs and the holes in a Mott insulator thermalize 
among themselves at a temperature T e higher than the 
lattice temperature Tl, by means of the electron-electron 
interaction (the interaction V defined in the previous sec- 
tion in the case of narrow-band or the Hubbard on-site 
U in the case of a wide-band Mott-insulator) . This in- 
teraction between the electronic degrees of freedom leads 
to processes such as impact ionization and Auger recom- 
bination, as well as carrier-carrier scattering. This equi- 
librium state is achieved as follows. 

The energy current of the absorbed solar photons is 
given by 

f°° e 3 
Jabsorbed = I^fe L d£ exp(^)-l (10) 

where f2 s is the solid angle under which the Sun is seen, 
and T is the Sun's surface temperature T ~ 5760 K. 
Here, the lower energy cutoff is set by the value of the 
Mott gap U*. 

First, there is the impact ionization process, and, then, 
the Auger recombination process, i.e., the inverse of the 
impact ionization process, shown as the reverse process 
of the one shown in Fig. [3] These two processes can be 
also denoted as follows 

DOi o DO; + D0 2 + H, (11) 

where DOi is the initial DO which creates the additional 
DO-H pair in the Mott-insulator denoted as DO2, H, by 
changing its energy- momentum state to become DO^. In 
addition, there is the analogous processes for high energy 
holes, i.e., 

Hr oHi+DO + H 2 . (12) 

Furthermore, there is the electron-electron scattering 
(due to the matrix element V defined in Section IIIII or 
the on-site Coulomb repulsion U), which leads to DO- 
DO, hole-hole, and DO-hole interaction processes: 

DO1+DO2 O DO'i+DO^, (13) 
H!+H 2 O H' x + U' 2 , (14) 
DO + H o- DO' + H'. (15) 

If these scattering processes are excluded, and we con- 
sider only impact ionization and Auger recombination, 
then, only those processes involving DOs and holes in cer- 
tain energy and momentum range which allow for impact 
ionization and its inverse through energy-momentum 
conservation would take place and everything else would 
be excluded. As a result, the final distribution of carri- 
ers would depend on energy and momentum and, thus, 
it cannot be described by a Fermi-Dirac distribution. 
Therefore, in the absence of electron-phonon interaction, 
the above scattering processes are important in order to 
establish a common temperature and a Fermi-Dirac dis- 
tribution. 
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If we now consider the equilibrium of all of the above 
processes, we realize that they produce a distribution of 
DO-H pairs in equilibrium with a distribution of photons 
over energy due to the emission process: 



DO + H 07. 



(16) 



in which a DO and a hole (H) combine and a pho- 
ton (7) is emitted (luminescence) or the inverse where 
a photon produces a DO-H pair in the Mott insulator. 
This distribution of photons in equilibrium with a gas 
of quasiparticles consisted of DOs and holes is analogous 
to the case of electrons and holes in a band-insulator 
which can be describe d 18 ' 29 by a thermal distribution of 
photons at a temperature T e an d a chemical potential 
M7 = Vdo-h = (see Refs.lHSl). The value of T e is 
determined from the equation 



J emitted — J absorbed) 

n 

J emitted 



4TT 3 h 3 C 2 



de 



' - M-y \ 



1 



(17) 
(18) 



Here Q e is the solid angle of the emitter and it is ir for a 
planar emitter. This equation defines a temperature T e of 
the electronic system which is considered decoupled from 
the lattice and the phonons. A value of udo-h = is 
expected for a system in which the particle number is 
not conserved as is the case for the DO and holes in 
the equilibrium state of impact ionization and Auger- 
recombination. 

Assuming that the density of the DO and of holes near 
the interface is not large, the above process produces an 
equilibrium Fermi-distribution of DO and holes at the 
above temperature T e with no separation of the quasi- 
Fermi-energies, i.e, the chemical potential difference be- 
tween DO and holes Udo-h is zero. 

Now, equilibrium is achieved (under open circuit con- 
ditions) at the temperature T e which is higher than the 
lattice temperature Tj,. As shown in the previous sec- 
tions, the time scale for the Mott insulator to reach equi- 
librium (through impact ionization, Auger recombination 
and electron-electron scattering) is much shorter that the 
time scale required for such a system to reach thermal 
equilibrium with the lattice through the process 



DO + H o r 



(19) 



where T denotes a phonon. This is so because of the facts 
that i) the time scale for the decay processes described by 
Eqs.fTTIand [T2lis much faster than the process of electron 
decay via phonons and ii) the interaction matrix elements 
involved in the processes described by Eq. [13l 1141151 are 
much stronger than the electron-phonon interaction ma- 
trix elements leading to the processes Q1J] These facts 
allow first equilibration of the electronic system to a tem- 
perature T e , which is different than Tl by creating mul- 
tiple DO-H carriers per incident photon. This has been 
argued in Refs. ITa - flTI to lead to a significant increase 
of the efficiency by means of carrier multiplicatio n 16 ! 17 



for a conventional solar cell. By treating the carriers in 
the Mott-insulator as weakly interacting quasiparticles 
which obey Fermi-Dirac statistics, which may allow us 
to apply the same assumptions and approximations used 
by Brendel et alJ^ to calculate the efficiency of a conven- 
tional solar cell, we also find that the maximum efficiency 
is bounded from above by the maximum of the following 
function of the voltage V and the Mott gap U* 



v(v,u*) = 


q v 9 


9(U*) = 






Ju* 


r(V,U*) = 


/'OC 




Ju* 


P — 

1 in — 


Jo 



g(U*) - £ r(V,W) 



p. 

1 in 

m(e)e 2 



exp(e/fcgT e ) - 1 ' 

m(e)e 2 
exp((e - qV)/k B T e ) - V 

^3 



"exp(e/fc s T s ) - 1 



(20) 
(21) 
(22) 
(23) 



In the above expression, apart from a multiplicative con- 
stant common in the numerator and denominator ex- 
pressions which cancels out, g(V) corresponds to the 
generation current, r(V,U*) to the recombination cur- 
rent and Pi n the total power carried by the sunlight; 
m(e) = min([e/U*],mo) is the number of DO-H pairs 
corresponding to each photon in dynamic equilibrium, 
where [e/U*\ stands for the integer part of e/U* and mo 
is a maximum allowed value of the carrier multiplication. 
Furthermore, q is the carrier charge and £ = n/Q s is the 
ratio of ir to the solid angle through which the Sun is 
seen from the Earth and the case £ = 1 corresponds to 
fully concentrated sunlight. If we follow Werner et alJ£ 
and Brendel et alJ^ and use a value of T e = 300K, we 
reproduce the same upper bounds for the solar-cell effi- 
ciencies as a function of U* as those reported in Fig. 1 
and Fig. 2 of Ref. [l6|. The calculated values for these up- 
per bounds for solar-cell efficiencies reach up to 85 % for 
fully concentrated sunlight and in the limit of vanishing 
value of U* when mo becomes very large. 



B. Further efficiency improvement 

As argued by Wiirfel™, however, since the phonons 
are excluded from the equilibration process, the DO and 
holes through Eq. [T7] will reach an equilibrium tempera- 
ture T e which should be higher than 300 K; as a result 
of the fact that the equilibrium temperature of the elec- 
tronic degrees of freedom is higher than 300 K assumed in 
the calculation of Brendel el alj^ 6 -, the efficiency should 
be lower than the values calculated using T e — 300 K. 
As was argued in the previous sections, the electronic 
equilibration is achieved much faster than the equilibra- 
tion with the lattice, and this is important for another 
reason, namely, it allows the application of an idea dis- 
cussed by Wiirfel^, which is summarized below, in order 
to improve the efficiency further. 

The photovoltaic cell can achieve maximum efficiency 
when there is a way to cool the gas of carriers to the 
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lattice temperature Tl via an isentropic processes where 
the relatively high energy of the carriers is not dissipated 
by heating the lattice but rather by creating electrical 
energy. Following Wiirfeli^ we assume that (a) impact 
ionization (and Auger recombination) occurs at much 
higher rate than photon absorption/emission, which as 
argued above is more likely in a Mott-insulator based 
solar-cell and (b) the carriers are removed from the solar- 
cell through membrane leads which are attached to the 
solar-cell on each side-iS and which can be semiconductors 
of n-type and p-type with narrow bands of band-width 
6e and large band-gaps. 

Because of the assumption (a) discussed above, the 
equilibrium of the carriers is not disturbed by the with- 
drawal of electrons and holes when the circuit is short. 
Now, the interaction with the slower acting phonons of 
the membranes cannot be avoided; in this case, however, 
this interaction is utilized in a favorable way. Namely, 
the hot electrons and holes with energies e e and in 
a narrow interval Se and with temperature T e can be 
cooled to the temperature of the lead-membranes Tl by 
keeping the entropy constant. This can occur due to 
the above mentioned assumption (b) where the carriers, 
whose concentration is constant, are transported via an 
approximately flat band, which prevents energy losses 
because the carriers keep approximately the same en- 
ergy. This works as a Carnot engine which produces work 
Heh = (1 — 7L/T e )(e e + eh) which appears as a difference 
of the quasi-Fermi levels between the electron and hole 
carriers in the two membranes and is given by qV. The 
details of how this works as well as the calculation of the 
efficiency can be found in Ref. [TH. This calculation leads 
to a maximum energy efficiency of 85 % also in the limit 
of a vanishing Mott-gap and for maximum concentration 
of solar cell radiation. Therefore, in this subsection, we 
conclude that the idea proposed by Wurfel may be uti- 
lized in solar-cells based on Mott-insulators. 



C. Other issues regarding efficiency 

A significant reason why a narrow-band Mott-insulator 
may be preferred over the wide-band case is that the 
threshold for impact ionization to occur for the latter 
case (and in a wide-band semiconducto r 1 ' 30 ! 31 ) is approx- 
imately 3 times the Mott-insulator gap. This is because 
energy and momentum are conserved together via the im- 
pact ionization (or Auger recombination) and presence of 
a wide-band implies a strong dependence of the energy 
with momentum. On the other hand, a narrow-band 
Mott insulator of the case discussed in Sec. [Hi] gives a 
much smaller energy threshold for impact ionization. 

Finally, there is an important issue to address which 
is that in a Mott insulator the carrier mobility might be 
low. However, the effective collection length, as well as 
the drift field of the minority carriers can be large in such 
a Mott insulator, due to the following reasons. The main 
reason why a DO diffuses from the n-doped region in- 



side the p-region to fill a hole, is the large energy gain 
which will be resulted by the fact that the double oc- 
cupancy in a Mott insulator of the type considered here, 
costs several eV. These materials become Mott insulators 
precisely because of this large on-site Coulomb repulsion 
which prevents the electrons from hopping from one site 
to a nearest neighboring site and this leads to an insulator 
in the half-filled band case. Therefore, this strong on-site 
Coulomb repulsion between electrons produces high dif- 
fusivity because the DO in the n-doped region move to 
the p-doped region to gain a large amount of energy by 
filling the hole and this may lead to a high drift voltage 
and large diffusion length regardless of the fact that the 
electron mobility might be low. 



VI. CONCLUSIONS 

What should be the nature of the absorption spectrum 
in the case of the Mott insulator considered here in or- 
der for it to be a candidate for the basis of an efficient 
solar cell? There are two cases as discussed in the paper 
for the case of a narrow gap Mott insulator, a) Narrow 
Mott-band materials. In this case in order for the effect 
considered in the present paper to be significant, we need 
other bands present, not necessarily of Mott type. The 
incident solar photon may be absorbed by an electron 
from another such valence band or the Mott conduction 
band and be promoted to a highly-excited band. If there 
are no issues related to selection rules, these electrons are 
expected to decay very quickly (for the reasons discussed 
in the paper) to produce one or multiple DO-hole pairs 
from the Mott band. Therefore, if the material is appro- 
priately selected, a large fraction of the solar spectrum, 
could be absorbed leading to a much larger photo-current 
than in band-semiconductors, b) In a wide-band Mott 
insulator, almost the entire solar spectrum could be ab- 
sorbed by the same single Mott-band. As analyzed in the 
present paper, in this case the initially photo-excited sin- 
gle DO-hole pair decays into two or more DO-hole pairs 
by means of the on-site Coulomb repulsion within a time- 
scale much shorter than other internal processes which 
dissipate the excess energy of the DO or hole beyond the 
Mott gap. 

In conclusion, it was argued that a p-n hetero junction 
in which the illuminated side is a narrow-gap Mott insula- 
tor, if fabricated, should be expected to behave as a high 
quantum efficiency solar cell. This can be achieved when 
the Mott insulator's energy spectrum for multiple quasi- 
electron-hole pairs falls within the energy range of solar 
photons. We find that in such relatively high- 

energy incident photon initially creates a high-energy 
electron-hole pair which quickly dissociates into multiple 
quasi-electron-hole pairs and this leads to an increased 
photo-current. As required for the device to work, it 
is estimated that this dissociation occurs in a time scale 
shorter than the time needed for the excited electron-hole 
pair to relax through the emission of phonons or other 
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intrinsic excitations. This allows equilibration of the elec- VII. ACKNOWLEDGMENTS 

tronic degrees of freedom among themselves, with excita- 
tion of multiple DO-H pairs for every absorbed photon, 
at a different temperature, on a time scale significantly 
shorter than the time required for the system to reach 
equilibrium with the lattice. This separation of time 

scales also allows the application of other previously pro- I wish to thank V. Dobrosavljevic, B. I. Halpcrin, C. S. 
posed ideas to improve solar-cell efficiency by cooling the Hellberg, P. Schlottmann, and S. von Molnar for useful 
carriers isentropically to a lower temperature^. discussions. 
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